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Abstract. The Heston stochastic volatility process is a degenerate diffusion process where the 
degeneracy in the diffusion coefficient is proportional to the square root of the distance to the 
(— I boundary of the half-plane. The generator of this process with killing, called the elliptic Heston 

^ operator, is a second-order, degenerate-elliptic partial differential operator, where the degeneracy 

in the operator symbol is proportional to the distance to the boundary of the half-plane. In 
mathematical finance, solutions to obstacle problem for the elliptic Heston operator correspond 
i | to value functions for perpetual American-style options on the underlying asset. With the aid 

Qh of weighted Sobolev spaces and weighted Holder spaces, we establish the optimal C 1 ' 1 regularity 

(up to the boundary of the half-plane) for solutions to obstacle problems for the elliptic Heston 
f"1 operator when the obstacle functions are sufficiently smooth. 
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1. Introduction 



In jl] , the authors established the existence and uniqueness of a solution in a weighted Sobolev 
spacd^ u € H 2 (G, to), to the obstacle problem, 

min{^4ti — /, u — = a.e. on G, 

(1.1) 



on r 



for the Heston operator [13J, 



Au:= -| (u^ + 2gau xy + cx 2 n ra ) - (r - g - |j ^ - «(0 - y)u y + rtt, (1.2) 

on a subdomain G (possibly unbounded) of the upper half-plane H := lx(0, oo), where / : G — > M. 
is a source function, g : T\ — > M prescribes a Dirichlet boundary condition along T\ :=Hfl dG ', 
and ip : G U Y\ — )• R is an obstacle function which is compatible with g in the sense that ip < g 



on T\. The differential operator A given in (1.2) is elliptic on G but becomes degenerate along 



To, where To denotes the interior of {y = 0} n dG. Because n0 > (see assumption ( |1.3| ) below), 
no boundary condition is prescribed along the portion To of the boundary (9^ = Tq U Ti of G. 

The operator A is the generator of the two-dimensional Heston stochastic volatility process with 
killing, a degenerate diffusion process well known in mathematical finance and a paradigm for a 



broad class of degenerate diffusion processes. The coefficients defining A in (1.2) are constants 
assumed throughout this article to obey 

o-^O, — 1< e < 1, r>0, q>0, k > 0, 9 > 0, (1.3) 

while their financial meaning is described in |13j . For a detailed introduction to the Heston 



operator and the obstacle problem (1.1), we refer the reader to our article [3]. 



In this article, we will establish d' 1 regularity on G U Tq and a priori Ci' 1 estimates for the 



solution u to (1.1) on subdomains U t£= GU Tq. We use C^' 1 to indicate a weighted Holder norm 
and corresponding Holder space which are distinct from the usual C 1 ' Holder norm and Holder 
space and which take into account the degeneracy of the operator, A, along y = — see section 
[2] for their definition. In the case of a uniformly elliptic operator on a bounded domain, interior 
C 1 ' regularity was established by Brezis and Kinderlehrer [2j (see also jl 1|, Theorem 1.4.1] for a 
statement of their result and an exposition of their proof), while global C 1 ' regularity, given a 
Dirichlet boundary condition, was established by Jensen [13J (see also \23\ Theorem 4.38] for a 
statement of his result and an exposition of his proof), recalling that [UJ p. 23], for a bounded 
domain U C W 1 , one has W 2,00 (U) = C l,1 (U). To the best of the authors' knowledge, however, 
our article is the first to establish C • regularity of a solution to an obstacle problem defined by a 
degenerate elliptic operator, despite the importance of this question in applications to American- 
style option pricing problems for asset prices modeled by stochastic volatility processes [13] . 

For interior C ' regularity, the case of a uniformly elliptic operator on a bounded domain 
reduces, by standard methods (see, for example, [U [23]), to the case of the Laplace operator 



^See section j2j for its definition. 
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and ingenious techniques introduced by Caffarelli [3] greatly simplify the proof of interior C ' 
regularity for solutions to an obstacle problem in this case. We shall adapt Caffarelli's approach 
in our article but, because of the degeneracy of our operator, A, along y = 0, careful consideration 
must be given to the different scaling of the equation near y = 0. This scaling is reflected in the 
use of the cycloidal distance function, s(-, •), defined in section [2] and of weighted Sobolev, Holder 
and C ' spaces. Weighted Sobolev and Holder spaces have been introduced previously (see, for 
example, [U [7J [TB]) in order to obtain sharp estimates for solutions to equations involving 
degenerate elliptic operators of the form (|1.2|) and their parabolic analogues. 



Let B p (Q ) := {Pel 2 : dist(P, Qo) < p} denote the open ball with center Qo = (po, qo) G M 
and radius p > 0, and set 

B+(Q ):=B p (Q )nU. (1.4) 

For a given radius Ro > and for any R obeying < R < Ro, we denote 

V:=B+ Q (Q ) and U := B+(Q ). (1.5) 

Throughout our article, we shall assume that Qo = (po,Qo) E S with < qo < A, for a positive 
constant]^] A. We shall abuse notation slightly and let To denote the interiors of dWndG, dWndV, 
or dWndU when we write GuTq, VUTo, or UUTo, respectively. The definitions of the weighted 
Holder spaces, Cf{V), C 2+a (V) and Cl'\V), which we require for the statement of the main 
result of this article below are collected in section [2j 

Theorem 1.1 (Optimal regularity). Let Ro > and A > and suppose Qo = (po,qo) £ S with 
< qo < A. Let V be as in (1.5). Assume that u E H 2 (V, ro) n C(V) is a solution to the obstacle 



problem, 

mm{Au — f, u — tp} = a.e. on V, (1-6) 

with ^ G C 2 s +a {V) and f G Cf{V), for some a G (0, 1). Then, u G C l s ' l {V U T ) and there is a 
constant C, depending on oc,Rq,A, and the coefficients of the operator A, such that if U is as in 



fll.5| ) with R = R /2, then 

\H\ c y(u) - c (lMlc(v) + ll/llcf(v) + IMIom(v)) • ( L7 ) 

Theorem |1.1| immediately yields 

Corollary 1.2 (Optimal regularity). Let & C H be a bounded domain. Assume that u £ 



H (G , ro) is a solution to the obstacle problem (1.6) on & with ip G C 2+a (G) and f G C, 



for some a G (0, 1). Then, u G C^' [ff U To) and, /or each precompact subdomain &' <g ^ U To, 
i/iere zs a constant C, depending on a, G 1 , G , and the coefficients of the operator A, such that 



Note that A in (1.2| is uniformly elliptic on Br (Qq) when go > -Ro, and results concerning regularity of 



solutions to ( 1.6 \ are then standard [111123] and so, for the purpose of this article, we could choose A = Ro without 



loss of generality. 
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Our proof of Theorem proceeds by adapting ideas of Caffarelli in [3| ; see also an exposition 
by Petrosyan in [19]. However, because our operator is degenerate, careful consideration must 
be given to the difference of the scaling of the equation in regions close (y small) and away 
(y > P > 0) from the portion of the boundary, {y = 0} n dV, where A becomes degenerate. 

1.1. Generalizations. When the main result of our article (Theorem |1.1[ ) is combined with 



Jensen's global C ' regularity theorem p3], we see that H (ff, to) solutions, u, to (1.1) actually 



belong to C s ' (0) n C ' {& U T\) under hypotheses on / and ip analogous to those stated in 
Theorem |l.l| By making further use of methods in Til] . it should follow that u S ^l'\ c(^' 
Moreover, there is good reason to believe that results on the regularity of the free boundary 
for the obstacle problem defined by a non-degenerate elliptic or parabolic operator extend to 
degenerate operators of the kind considered in this article; see |20j and references therein for the 
non-degenerate elliptic case and \17\ fl"8] and references therein for the non-degenerate parabolic 
case. We shall leave consideration of these extensions to our future articles. 



The solution, u, to ( 1.1 ) can be interpreted as the value function for a perpetual American-style 
option with payoff function, ip [15]. The C^ +a {V) regularity property assumed for the obstacle 
function, ip, in the statement of Theorem 1 1 . 1 1 does not reflect the more typical Lipschitz regularity 
for if) encountered in applications to mathematical finance, such as tfj(x,y) = m&x{E — e x ,0}, 
where E is a positive constant, in the case of a put option. Nevertheless, simple examples in this 
context [221 §8.3] and results of [TT] [18] suggest that the solution, u, should nevertheless have the 
optimal Cs' 1 regularity even when ip = max{E — e x ,0}. Again, we shall leave consideration of 
this question to our future articles. 

We have chosen, in this article, to work with our model, the Heston operator A, because of 
its relevance to mathematical finance and reliance on results in our previous work [1] and that 
of Feehan and Pop [U [7J QUI EI] . However, we expect that the Ci' 1 regularity result and a priori 
estimate in Theorem 1 1 . 1 1 may be easily generalized to higher dimensions and degenerate elliptic 
operators on W 1-1 x (0,oo) with variable coefficients, 

under the assumptions that (a^-) is strictly elliptic, b n > v > 0, for some constant v > 0, and 
c > and all coefficients are Holder continuous of class C"(V), for some a 6 (0,1). See [8j for 
an analysis with applications to probability theory based on parabolic operators of this type. 

1.2. Outline of the article. For the convenience of the reader, we provide a brief outline of the 
article. We begin in ^2] by reviewing our definitions of weighted Holder spaces [5] and weighted 
Sobolev spaces [3] which we shall need for this article. In we review results from [6] [TJ [TO] 
concerning existence, uniqueness, and regularity of solutions to the elliptic Heston equation on 
bounded subdomains of the upper half-plane; see also [5]. In Q we develop the key pointwise 



growth estimates (see Propositions 4.1 and 4.4) for solutions to the obstacle problem for the 



elliptic Heston operator. We conclude in §[5] with the proof of our main result, Theorem |1.1[ 
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1.3. Notation. Throughout the rest of the article we will set Lu := —Au, where A is given by 



(1.2) and we work with L instead to facilitate comparisons with the methods of Caffarelli [3j and 



the sign conventions therein. The operator L is then given by 

Lu = | (u xx + 2gau xy + a 2 u yy ) + (r - q - |j u x + k (9 - y) u y - ru, (1.9) 



with coefficients which satisfy the assumption ( |1.3[ ). 

We let C = C(*, ...,*) denote a constant which depends at most on the quantities appearing 
on the parentheses. In a given context, constants denoted by C,C',--- and so on may have 
different values depending on the same set of arguments and may increase from one inequality 
to the next. Constants with values denoted by K, K' , • • • and so on are reserved for quantities 
which remain fixed. We let C(L) denote a constant which may depend on one or more of the 
constant coefficients of the operator L (that is, q, r, k, 6, g, a). 

1.4. Acknowledgments. We are grateful to Arshak Petroysan for sharing Mathematica code 
from his lecture notes |19j and which we adapted to create the figures in this article. We are also 
grateful to Camelia Pop for many helpful conversations. 

2. Weighted Sobolev and Holder spaces 

In [2] the authors defined the following weighted Sobolev spaces of functions on a possibly 
unbounded domain ^Cl. 

Definition 2.1 (Weighted Sobolev spaces). Let L 2 (G, ro) denote the Hilbert space of Borel 
measurable functions, u : G — )■ K, such that 

1/2 



MlL 2 (^,re) : = (^j u 2 Vodxdy^j < oo, 



with weight function ro(x, y) := y@ 1 e 7 ' a ' for (x, y) G H, where /3 := 2k6 /a 2 and fi := 2k/g 2 
and the constant 7 > depends only on the coefficients of A. We define the vector space, 

H 2 (0,xo) := j-u G L 2 (tff,to) : y\D 2 u\, {1 + y)\Du\, (1 + y) l/2 u 6 L 2 {0, ro)} , 

where Du = {u x ,u y ) and D 2 u = (u xx ,u xy ,u yx ,u yy ) are defined in the sense of distributions. 

When equipped with the norm, 

1/2 



\ U \\H 2 (&,K>) : = 



J {y 2 \D 2 u\ 2 + {l + y) 2 \Du\ 2 + {l + y)u 2 ) todxdy^j , 



one finds that H 2 (i?,)x>) is a Hilbert space and, noting that G has dimension two, H 2 (G,to) C 
C(GuTi) via the embedding theorem for standard, unweighted Sobolev spaces [I], but elementary 
examples show that functions in H 2 (G, tn) need not be continuous up to f o- 
We next define weig hted C 1 ' 1 and Holder norms on a bounded domain ^Ci. 
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Definition 2.2 (Ci' 1 norm and Banach space). We say that u G Cl' l (&) if u belongs to C 1,l {ff)r\ 
C x (d) and 

:= l|y^ 2 «IU°°(^) + \\ D u\\c(ff) + \\ u \\c(0) < 00 • 

Also, we say that u G C s ' (ff U To), if u G Cs' (?7) for any subdomain U <s ^ U IV 

We recall the definition of the distance function, s(-, •) on H, equivalent to the distance function 
defined by the cycloidal metric, y~ 1 (dx 2 + dy 2 ) on M, and introduced by Daskalopoulos and 
Hamilton in [5] and by H. Koch in |16j . 

s(z,z ):=— Vz = (x, y),ZQ = (xq, yo) G H. (2.1) 

y/y + Vm> + v l x ~ x o\ + \y- vol 



This is the natural metric for our degenerate equation; see [5] for a discussion. The following 
weighted Holder spaces were introduced by Daskalopoulos and Hamilton in [5]. 

Definition 2.3 (C? and C 2+a norms and Banach spaces). Given a G (0,1), we say that u G 
C?(ff) if u G C(<?) and 



|u(zi) - u(z 2 )| 



Mlc?(<?) : ~ IMIc(^) + SU P s ( Zl < °°- 



We say that u G C 2+Q (^) if it has continuous first and second derivatives, Du,D 2 u, in and 
-Dn, yD 2 u extend continuously up to the boundary, dff, and the extensions belong to Cf(0). We 
denote 

ll n llc? +a (^) := ll n llc?(«?) + ll-^MIc? (d) + \\yD 2u \\cf {dy 

We say that|^]u G Cf(0 U Tq) if u G Cf(U) for every subdomain [/ <s U To and similarly that 
u G C s 2+a (^ U T ) if n G C s 2+a ([7) for every subdomain U m & U T . 



One can show that C 2+a (&) is a Banach space when equipped with the norm in Definition 2.3 

Remark 2.4. On any bounded subdomain [/ C EI we have, 

c\z- z Q \ < s(z,z ) < y/\z - z \, (2.2) 

for some positive constant c := c(diam(?7)) depending only on the diameter of U. Hence, C a (U) C 
Cf(U) C C a ' 2 {U). 



In Op. 901], when defining the spaces C?(sf) and C'i; +a {£/), it is assumed that si is a compact subset of the 
closed half-plane, {y > 0}. 
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3. SCHAUDER EXISTENCE, UNIQUENESS, AND REGULARITY RESULTS 

We collect some known results for solutions to the degenerate elliptic equation, 

Lv = f on V, (3.1) 
where V is as in (1.5). These results will be used in the proof of Theorem |1.1| Theorems 3.1, 3.2[ 



3.3 and |3.4| are proved in [10] and may be viewed as analogues of [TH Theorems 6.2, 6.6, 6.13 



and 6.14] and a generalization of [5] Theorem 1. 1.1]. 

Theorem 3.1 (A priori Schauder interior estimate). [10] Given f £ Cf{V), where V is as in 
flL5l ) ; and a solution u £ C 2 s +a (V U T ) n C(V) to 

Lv = f on V, 

there is a constant, C, depending at most on a, R, Rq, A, and the coefficients of L, such that ifU 
is as in (1.5), then 



v\\ C 2+* { u) <C(\\u 



\C(V) 



+ 



Cf{V) 



(3.2) 



Theorem 3.2 (Existence of a solution to a Dirichlet problem with continuous boundary data). 



[TO] Given f £ Cf{V UTo), where V is as in (1.5), and g 6 C(WndV), there exists a unique 
solution u £ C 2+a (V U T ) Pi C(V) to the Dirichlet problem 



Lv = f on V and v = g on M D dV. 



(3.3) 



Theorem 3.3 (A priori global Schauder estimate). [JU] Given f £ C?(V) and g £ C 2 s +a {V) and 
a solution u £ C 2+a (V) to the Dirichlet problem (3.3), there is a constant, C , depending at most 
on a,Ro,A, and the coefficients of L, such that 



Mlc s 2+Q (u) - c (ll n llc(u) + ll/llc? (v) + \\9 Ilc2 +a (v)) 



(3.4) 



Theorem 3.4 (Existence of a solution to a Dirichlet problem). [10] Given f £ Cf(V) and 
g £ C 2+a {V), there exists a unique solution u £ C 2+a (V) to the Dirichlet problem (3.3). 

The preceding results easily imply the following consequence when combined with a regularity 
theorem from [9] and a maximum principle estimate from [6]. 

Proposition 3.5 (Regularity and interior Schauder estimate). Let f £ C"(V U To) and let 
v £ H 2 (V, to) be a solution to 

Lv = f a.e. on V. 

Then, v £ C 2+a (V UTq). Moreover, if f £ Cf(V), there is a constant, C, depending at most on 
a, R, Rq, A, and the coefficients of L, such that if U is as in (1.5), then 

+ 



\v\\ C 2+ a(0) < c ( |M 



C(V) 



Cf(V) 



(3.5) 
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Proof. Choose R 1 obeying R<R X <R and let Vi := (Q ), so that U £ Vi <e V U T . Then 
/ G (Vi) and we may choose w G Cg +Q (Vi) to be the unique solution to Lw = / on V~i and 
w = on EI n <9Vi provided by Theorem 3.4 Moreover, vq := v — w £ H 2 (Vi, to) is a solution to 
Lu = a.e. on V\ and so, by [9j, we have v G C°°(ViUro) and thus v = v + w G C 2 s +a {Vi UT ). 
Since i?i is arbitrary, we obtain u G Cg +a! (V U Tq), as desired. Now fix R\ = {R + i?o)/2 and 



observe that (3.2) gives 



\ v o\\c* +a (U) - C \\M\c{v^)- 



3.4 



implies 



Setting <8i := {u G C£ +a (Vi) :ti = 0oni(1 dVy} and «8i := then Theorem 

that L : 5Si — )• 5^2 is bijective and bounded and thus L^ 1 : 5^2 - ^ ©l is bounded by the open 
mapping theorem. Hence, there is a constant, C, depending at most on a,R, Ro,A, and the 
coefficients of L, such that 



\w\ 



c s 2+a (v\) - C W L 1 



Cll/ll 



C»(Vi)- 



Therefore, by combining the preceding estimates, 



Mcg+uft < \\v \\ c 2+ a{U) + \\w\\ c 2 +a{0) 



<C[\\v 



\C(Vi) 



+ 



< C f llw||c(v\) + Ikllr-n?^ + 



< C \\v 



IC(Vi) 



+ 



C-(Vi) 
b(Vi) + \\J HCf(Vi) 



and thus (3.5) follows. 



□ 



Remark 3.6 (Alternative proofs of regularity in Proposition 3.5). The regularity assertion, v G 



Cg Jra (VUTo), also follows by using Theorem 3.2 instead of Theorem 3.4 in the proof of Proposition 



We can avoid relying on the regularity results in [9] if we are given v G H 2 (V, to) n C(V). 

provides a unique solution v G C 2+a {V U To) n C(V') to = / on V and 
n dV. But then u G -ff 2 (V, tt)) and by the weak maximum principle for L acting on 
functions in H 2 (V, to) [H Lemma 6.13 & Theorem 8.8], we must have v = v a.e. on V and thus 



Indeed, Theorem 3.2 
v = v on 



v G Ci +a (VUT Q )nC(V). 



Remark 3.7 (Alternative proof of the interior Schauder a priori estimate in Proposition 3.5). We 
can instead choose w G C 2+a (V U Tq) n C(V") to be the unique solution to Lw = f on V and 
on in dV provided by Theorem 3.2 If we also assume that the coefficient r in (1.9) is 



w 



positive, then the weak maximum principle estimate [61 Proposition 2.19 & Theorem 5.1] gives 

1 
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C(V) 



while (3.4) yields 



\w\ 



C 2 + a {U) 



< C 



< 



w 



c{vy 



C(V) 



+ 



Cf(V) 
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We now obtain v E C 2+a (V U Tq) n C(V) and the interior Schauder a priori estimate (3.5) by 
combining the preceding observations for uu with those for vq in the proof of Proposition 3.5 but 



without a need to appeal to Theorem 3.4 



The following weak and strong maximum principles are shown in [6j. Recall that if v E 
C 2+a (GU r ), then Dv E C(GU Tq) and yD 2 v E C(GU Tq) (by definition) while yD 2 v = on 
Tq (see [5j Proposition 1.12.1] or [HJ Lemma 3.1]). 

Theorem 3.8 (Weak maximum principle for the Heston operator). [6j Theorem 5.1] Let v E 
C 2+a (GUTo) flC(^) be a subsolution, Lv > on G and v < on WndG, for a bounded domain 
G C H. Then v < on G. 

Theorem 3.9 (Strong maximum principle for the Heston operator). [61 Theorem 4.3] Let v E 
Cl +a {ffUTQ)C\C{0) be a subsolution, Lv > on G, for a bounded, connected domain ^Ci. If 
v achieves its maximum value at a point P E G U Tq and, in addition, v(P) > if r > (where 
r is the coefficient of L in (1.9),), then v must be a constant on G. 

We finish this section by showing how to reduce to the case / = in Theorem |1.1| 

Proposition 3.10 (Reduction to a homogeneous obstacle problem). We may assume, without 
loss of generality, that f = on V in Theorem \1.1\ 

Proof. Let v E C 2+a (V) be the solution to the Dirichlet problem Lv = f on V and v = on 
in dV (its existence follows from Theorem 3.4). It follows from (3.4) that 

IMIc2+ a (y) - c (iMlc(v') + ll/llc<*(y)J ) 



and hence (see the proof of (3.5)) 



Mlc?+ a (V) - C \\f\\c*{V)' 



(3.6) 



where in (3.6) we use C to denote a constant which depends at most on a, Rq, A and the 
coefficients of L. 

Now if u is a solution to the obstacle problem (1.1) on V as in Theorem |1.1[ then 



u := u — v 



is a solution to the obstacle problem (1.1) on V with source function / = on V and obstacle 
ij) := ip — v on V. If Theorem 



1.1 



is proved for / = in place of / on V, then u E C s ' (U) and 



the estimate (1.7) for u yields 



\ u \\cy(u) ^ c 



lew 



+ 



But u = u + v E Cl' l (U) and we obtain the estimate (1.7) for u from the preceding inequality 
and the estimate (3.6) for v. □ 



Because of the reduction in Proposition 3.10 we may assume without loss of generality that 
u E H 2 (V, to) n C(V) is a solution to the obstacle problem (1.1) with obstacle function ip E 
C 2+a (V) and / = on V, that is, 

min{Ln, u — tp} = a.e. on V. (3-7) 
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We make this assumption for the remainder of this article. 



4. SUPREMUM BOUNDS 



We will assume, throughout this section, that u is a solution to the obstacle problem (3.7) 



on V, where V is as in (1.5), and that all the assumptions of Theorem hold. Adopting the 
terminology of mathematical finance, we call 



= {P G V U T : u(P) > ^b(P)} 
the continuation region (or non-coincidence set), 

g{u) = {P G V U T : u(P) = ip(P)} 
the exercise region (or coincidence set), and 

j?( n ) = (y u r ) n <9^(u) 



(4.1) 



(4.2) 



(4.3) 



the free boundary (or optimal exercise boundary, as it is known in mathematical finance). From 



(3.7) and (4.1), we see that 



Lu < a.e. on V and Lu = on ^(u) 



(4.4) 



3.5 



that u is of class C^ +a on ^(u). (Actually 



Since Lu = on ^(u), it follows from Proposition 
one may also easily see that u is of class C°° on ^(u).) 

We will establish sharp growth estimates from above on u — ip near free boundary points 
Pq G j^(it). Because of the degeneracy of our operator L, we will need to scale our estimates in 
different ways, depending on the distance of Po from the boundary portion, & n dW = {y = 0}. 
Similar estimates in the non-degenerate case, where L is the Laplace operator, A, were established 
by Caffarelli in [3]. 

The first such estimate, in Proposition 4.1 concerns with free boundary points Pq = (xq, yo) G 



&{u) with yo > 0. To simplify the notation we will assume that < yo < 1- The estimate near 
any free boundary point Po = (xq, yo) G J^"(u) with yo > 1 can be shown similarly. We have the 
following analogue of (3J Lemma 2]; see also [191 Lemma 1.6] (where L = A and ip = 0). 

Proposition 4.1 (Quadratic growth of solution near free boundary and away from degenerate 
boundary). Let u be as in Theorem \l.l\ and let Pq = (xo,yo) G &(u) D V with < yo < 1. Then 
there are constants < po < 1 and < C < oo, depending at most on the coefficients of L, such 
that if B poyo (Po) (e V , then 



sup (u - ip) < Cyop" 

Z PV0 /2(P0) 



Vp < po- 



(4.5) 



Remark 4.2. We shall establish (4.5) with the aid of certain auxiliary functions, £ in (4.9) and w 
in (4.19), defined on balls B pyo (Po). 
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Figure 4.1. Regions in the proof of Proposition 4.1 for estimating the growth of 
a solution near the free boundary and away from the degenerate boundary. 

We begin by observing that since B poyo (Po) (g V by assumption, the operator L is uniformly 
elliptic on B poyo (Po) <g H. Consider the linear approximation, 

l Po (x,y) :=i;(P )+Dij(P )-(x-x ,y-y ), (i,y)£R 2 , (4.6) 

to our obstacle function ip at Pq. A direct calculation shows that 

\L(l Po )\<M onB poyo (P ), (4.7) 



where, noting that < po < 1 and < yo < 1 as in the hypotheses of Proposition 4.1 and that 

< 



CHB po y {P )) 



M := K 



C M (*«>»o(flO)' 



(4.8) 



and the constant K > depends at most on the coefficients of L. For < p < po, let £ 6 
C 2+Q (-B p j /0 (Po)) be the unique solution (assured by [HI Theorem 6.14]) to the elliptic boundary 
value problem, 

L( = L(l Po ) on Bpy (P ), 
( = WMy Q p 2 ondB pyo (P ). 
The next lemma provides sharp bounds from above and below on £ in terms of p and the constant 



(4.9) 



M in (4.9). 



Lemma 4.3 (Quadratic growth of an auxiliary function near free boundary and away from 
degenerate boundary). The function ( G C 2+c "(B pyo (Po)) in ( |4.9[ ) satisfies the bound, 

My p 2 < ( < UMy p 2 on B m (P ), < p < p , (4.10) 



where M is as in (4.8) and po < 1 is a constant depending at most on the coefficients of L. 
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Before proceeding to the proof of Lemma 4.3 , we consider the effect of rescaling on the operator 
L. Observe that, for any u 6 C 2 (H), if 



v(x, y) := u(x + y x, y + y y), (x, y) := (x + y x, y + y y), 



then 



(Lu)(x,y) 



and therefore, 



yo + yoy _ 2 



y {vxx + 2gav xy + a v yy ) (x,y) + ( r - q 
+ k(9- yo(l + y)) j/) - rw(x, y), 



y (Lu)(x, y) = (L yo v)(x, y), V(x, y) G 



yo(i + y) 



2/0 V(»,y) 



where 



(4.11) 



(4.12) 



+ k (6 - y (l + y)) y) - ry v(x, y), V(x, y) € E. 
We now proceed to the 

Proof of Lemma\4.3\ Since the ellipticity constant for L depends on yo, we shall use the rescaling 



(4.13) 



in (4.11). Note that the operator L yo is uniformly elliptic on Bi/ 2 , since 

4 < ~2~ < 4 ° n ^Z 2 ' 
and the coefficients of L yo are bounded by a constant (recall that yo < 1) depending at most on 
the coefficients of L. Let 



C{x, y) ■= -^—C,{xq + y x, y + y y), V(x, y) £ B p , 
My 



with £ as in (4.9). It follows from (4.7) and (4.9) that £ satisfies 

\L yo C\ < 1 on Bp, 
C = 10p 2 on dB 



(4.14) 



(4.15) 



since (4.11 ) yields 



We will show that 



(L yo ()(x,y) 



j^m(LO^y) = ^(LC)(x,y). 



P 2 < C < 14/> 2 on Bp, 



(4.16) 



provided p < po, with po < 1 a constant depending at most on the coefficients of L, and this will 
conclude the proof of the lemma, since (4.10) follows from (4.14) and (4.16). 



To this end, we consider the barrier function, 

•&{x,y) := ax 2 , V(x,y)eIR 2 , 



(4.17) 
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for different choices of constants a£l and compute that 

2/o(l + y) 



(L yo $)(x,y) = a 



(1 + 1/) + 2 r 



x - ry x 



Since 1/2 < 1 + y < 3/2 on By 2 , by choosing p < po, with po < 1 a constant depending at most 
on the coefficients r,q of L, and using (x,y) £ -B p and recalling that < yo < 1) we can ensure 
that 

Ly Q -d < - a on B p , if a < and > - a on if a > 0. (4-18) 



Choose a = —8 and set w := Q + ~d — p 2 . By combining (4.15) and (4.18) and using the definition 
(4.12) of L yo , we obtain 

L yo w < 1 — 2 + rp 2 yo < on B p , 



if p < 1/r (remember that yo < !)■ On the other hand, since C = 10 p 2 on dB p by fl415| ) and 
using (4.17), we see that 

to = £ + t? - p 2 

> 10/j 2 - 8p 2 - p 2 > on 

Therefore, the weak maximum principle for L yo on B p implies that 

( + tf-p 2 >0 on B p . 

Since •& = —8x 2 < on I 2 , we conclude that ( > p 2 — 1!) > p 2 in B p . 

We will now estimate Q from above. This time we take a = 4 and setting z := £ + #, we now 
find from ( |4.15| ) and ( |4.18D that 

Lj,,,*: > -1 + ^4 = on 5 p . 



But (4.15) and the definition (4.17) give 1? = 4x 2 < Ap 2 on i? p and 

z = ( + ^ < 10 P 2 + 4p 2 = 14p 2 on dB i 



p- 



and so the weak maximum principle for L yo on B p shows that Q < 14p 2 on B p . This finishes the 
proof of (4.16), and hence concludes the proof of our lemma. □ 



Proof of Proposition J^.l. We shall follow the proof of Lemma 2 in [3]. Our case is more difficult 
since linear functions are not solutions to the equation Lu = 0. In addition, our operator L has 
variable coefficients and our scaling depends on the ellipticity constant of the operator L on V, 
which is comparable to yo- 

With l Po given by (|46j) and < p < p and C G C 2+a (B pyo (P )) the function defined by ( [4~9] ), 
we set 

w := u - l Po + C G tf 2 (5 rao (P )) n C(B pyo (P )) (4.19) 

and observe that 



w 



(u-iP) + {iP-l Po + ()>iP- l Po + C on Bpy^Po), 
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since u > ip on B pyo (Po). By Taylor's theorem, 



x,y) -l Po (x,y)\ < 2yop 2 U\\c^(B pyo (P ))i V(x,y) G B pyo (P ). 



(4.20) 



Since £ > Myop by (4.10) and yo < 1 (by hypothesis in Proposition 4.1), we conclude that 

„2 



w > -2 y y 

> -2y 2 p 2 



J P0V0 \ 



+ My p z 
+ My p 2 (since p < p ) 



> on Bpy Q (P ) (by definition of M), 



provided that the constant K in the definition (4.8) of M is chosen large enough that K > 2. 
Also, since L( = L(lp Q ) by (4.9), we have 

Lw = Lu < a.e. on B pyo (PQ), 



where the inequality follows from (4.4). Let us now split w as 

W = W\ + W2 , 



(4.21) 



where w\ G C 2+a (B p y Q (PQ))r\C(B pyo (Po)) is the unique solution (assured by [12j Theorem 6.13]) 
to 

Lwi =0 on B pyo (P ), 
wi = w on dB pyo (P ). 
(Note that w = u — lp + C belongs to C(dB pyo (Po)).) By the weak maximum principle, we have 

< wi < w on B pyo (P ), 

(since w > on B pyo (Po)), and hence 

< W2 < w on B pyo (Po). 



(4.22) 

e have 
(4.23) 

(4.24) 



The inequality (4.23) obeyed by w\ and the definition (4.19) of w yield, 

MPo) < w(P ) = ((Po), 

Wl {P ) < UMy p 2 . 



and thus, by (4.10), 



Consider the rescaled solution, 

wi(x,y) := w^xo + yox^o + yoy), V(x,y) G B p , 



(4.25) 



(4.26) 



and observe that the function w\ G C 2+a (B p )nC(B p ), by (4.11) and ( |4.22 ), satisfies the uniformly 
elliptic equation, 

L yo w\ = on B p . 



The Harnack inequality |121 Corollary 9.25 & Equation (9.47)], the definition (4.26) of W\, and 
the inequality (4.25) imply the estimate, 

sup m < C inf m < C"u)i(0) = C' Wl (P ) < CMy p 2 , 



B 



p/2 



B 



p/2 
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for constants C and C = 14C" which depend at most on the coefficients of L, but are independent 
of yo, and the constant M is given by (4.8). Hence, by (4.26), 

sup wi < CMy Q p 2 , < p < p . (4.27) 

B py /2( P °) 

We will next bound W2 on B pyo (PQ), taking care to note that (like the regularity of u in Theorem 
1.1) W2 only belongs to H 2 (B pyo (P )) n C(B pyo (P )). Recall that < W2 < w on B pyo (P ) by 
( |4.24[ ) and that w 2 = on dB pyo (P ) by ( |4.21[ ) and ( |4.22| ). Assume that Pi = (xi,?/i) is a 
maximum point for the function u>2 on the closure of the ball B pyo (Po) and that W2(Pi) > 0. 
Then, Pi G P rao (Po) and we consider two cases. 

Case 1 (Pi G ^(w))- If Pi G (where u = V), then u(Pi) = ?/>(Pi) and hence, by the 

inequalities ( |4.10[ ), ( |4.20[ ), ( |4.24[ ), and definition ( |4.19[ ) of w, we have 

MPi) < w{Pi) = HPi) - lp (Pl) + ((Pi) < 16My p 2 , 



provided the constant K in the definition (4.8) of M is chosen large enough that K > 2. 

Case 2 (Pi G ^(u)). If Pi G ^(n) (where u > ip) then, since Lu>2 = on the open set 
^(u) n B PVo (Pq) and u>2 achieves an interior maximum there, the strong maximum principle |12| 
Theorem 3.5] implies that W2 must be constant on the connected component of ^ (it) n B pyo (Po) 
containing Pi. Since W2 = on dB pyo (Po) and w 2 (Pi) > by assumption, it follows that 
^2 (Pi) = w 2 (P2) for some point Pj G £{u) n B /yyo (Po). (Recall that, by hypothesis, Pq G &(u) 
and so S{u) n B pyo (P ) is non-empty.) Thus, by the inequalities ( |4.10[ ), ( |4.20[ ), ( |4.24[ ), and 
definition (4.19) of u>, we have 

^ 2 (Pl) = ^ 2 (P 2 ) < ^(P 2 ) = V ) (^ 2 ) - lPo(P2) + C(P2) < 16Myop 2 , 



provided the constant K in the definition (4.8) of M is chosen large enough that K > 2. 
By combining the two cases and recalling that W2 < w 2 (Pi) on B pyo (Po), we obtain 



sup u>2 < 16Myop 

B pyo / 2 (Pb) 



0< p< po. 



(4.28) 



By combining the supremum bounds (4.27) and (4.28) for w\ and IU2, respectively, we obtain 

(4.29) 



IV 



<CMy p z on B pyo/2 (P ), < p < p , 



where C depends at most on the coefficients of L, and M is given by (4.8). This shows, in 
particular, by ( |4.10[ ) and ( |4.19[ ), that 

u-Ip < CMy p 2 on P p2/o/2 (Po), 



where C depends at most on the coefficients of L. Now, again using (4.20), we have 

u-tp = u- l Po + l Po < CMy p 2 on B m / 2 (P ), < p < p , 

for a possibly larger constant C that depends at most on the coefficients of L, and this gives the 
desired bound (4.5). □ 
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We will next establish a supremum bound for the solution, u, which holds near y = and is 
independent of the yo coordinate of the point Po- 

Proposition 4.4 (Linear growth of solution near free and degenerate boundaries). Let u be as 
in Theorem \l . 1\ and let Po = (xo,yo) G ^(u) n V with < yo < 9/4, where 9 > is a coefficient 
of L in ( |1.9[ ), Then, there are a constant < po < 1 and a constant < C < oo, depending at 
most on the coefficients of L, such that if B~^ q {Pq) <s yuTo, then 

sup (u-i/>) < Cp|MlcrM(B+ (Po))' 0<p<p - (4.30) 

B p + / 2 ^o) 




Figure 4.2. Regions in the proof of Proposition 4.4 for estimating the growth of 
a solution near the free boundary and near the degenerate boundary. 



Our proof of Proposition 4.4 follows the pattern of the proof of Proposition |4.1| However, we 
shall use a different scaling. Observe that 



\I4\<Nk$ on B+(P ) 



where k > 0,9 > are coefficients of L in (1.9) and 



N := K' 



Cl,l(B+(P ))' 



(4.31) 



(4.32) 



where K' is a constant which depends at most on the coefficients of L in (1.9) (remember that 
< y < 9/4 and < p < 1). 

For < p < po, let £ G C^ +a (B+(P )) D C(B+(P )) be the solution to the boundary value 
problem, 

'ht = L^ on B+(P ), 
£ = WNp onUndB p (Po), 



(4.33) 



provided by Theorem 3.2 or 3.4 
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Lemma 4.5 (Linear growth of an auxiliary function near free and degenerate boundaries). The 
function £ given by (4.33) satisfies the bound, 

Np<£<20Np on B+(P ), < p < p , 



(4.34) 



where N is as in (4.31) and po < 1 is a constant depending at most on the coefficients of L. 
Proof. We first establish the bound from above. We set 

z := £ + 2N(y - y + p) - 20iVp G C* +Q (B+(P )) D C(B+(P )), 



and use (1.9) to compute that 

Lz = L£ + 2Nk(6 -y) - 2Nr(y - y + p) + 20Nrp 



> -Nk6 + 2Nk{9 - y) - 2Nr(y - y + p) + 20Nrp (by fl43l| ) and ( |433| » 

> -JVk;0 + iVK6» + 16iVrp 
>0 on B+(P ), 

if < yo < 0/4 and p < po with po < min{6>/4, 1} and noting that 0<y<yo + P< #/2. On the 
other hand, since £ = lOiVp on Bin dB^(Po) by (4.33), we have 

z = £ + 2iV(y - y + p) - 20Np 
< WNp + 4Np - 20Np 
<0 onIflSB+(P ). 



Hence, the weak maximum principle for L on B^(Pq) (Theorem 3.8), implies that z < on 
B+(Pq), which implies the desired upper bound in (4.34), 

i = z- 2N(y -y + p) + 20Np < 20Np on B+(P ), 

since y-y + p>-p + p = 0on B+(P ). 
For the bound from below, we now set 

z := £ - 4N(y - y + p) - Np, 



and use (1.9) to compute that 



Lz = Li- 4Nk{6 -y) + ANr{y - y Q + p) + Nrp 



< NkO - 4Nk(9 -y) + ANr{y - y + p) + Nrp (by | |43i| ) and ( |433| |) 

< 7Vk(9 - 2JV/cfl + 9iVrp 

< 7Vk6I - 2Nk8 + A^k6» 
<0 onB+(P ), 
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if < yo < 0/4 and p < p with p < min{6>/4, Kd/(9r), 1} and noting that 0<y<y + p< 9/2. 
On the other hand, since £ = 10iVp on Bin <9P+(Po) by (4.33), we have 

z = i- AN(y -y + p)-Np 
> WNp - 8Np - Np 
>0 on Hn<9P+(P ). 



The weak maximum principle for L on P+(Po) (Theorem 3.8) once more shows that z > on 
P+(Po). We conclude that 

£ = z + 4N(y-y + p) + Np> Np on B+(P ), 

provided that p < po, with po < 1 depending at most on the coefficients of L. This yields the 
desired upper bound in (4.34) and finishes the proof of the lemma. □ 



We will now give the proof of Proposition 4.4 



Proof of Proposition 4-4- We give an argument which is similar to the one used in the proof of 



Proposition 4.1 but we scale our estimate differently and use Lemma 4.5 instead of Lemma |4.3| 
We set 

»:=«-^ + (£ P 2 (P+(P ), to) n C(P+(P )), (4.35) 



(4.36) 



with £ G C* +a (B+(P )) n C(B+(P )) given by ([433]). Then w satisfies 

Lw = Lu a.e. on P+(Po). 

Let us now split w as w = w\ + w 2 , where w\ G C 2+a (P+ (Po)) H C(P+(Po)) (whose existence is 
assured by Theorem 3.2) is defined by 

Lwi = on B+(P ), 



w x = w on Mn dB+ (P ). 



(4.37) 



By the weak maximum principle for L on P+(Po) (Theorem 3.8) and the fact that w > £ > on 
P+(Po), we have 



< w 1 < w on BJ (P ), 

and thus 

0<w 2 <w on B+(P ). 
From ( j4.34| ), ( |4.35| ), ( |4.38| ), and the fact that u(P ) = ip(Po), we see that 

Wi(P ) < w(P ) = £(P ) < 20Np on S+(P ). 
Set (x p , y p ) := (xo + py, yo + py) and consider the rescaled solution, 

wi(x,y) := w 1 (x + px,y + py), (x,y) G Pi n {y p > 0}, 
which satisfies the equation, 

L p w 1 = onPin{y p >0}, 



(4.38) 
(4.39) 

(4.40) 



(4.41) 
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where (compare (4.12)) 



— [v. 
+ k( 



+ 2gav xy + cr 2 Vyy) + (r-q- y p )v x 



(4.42) 



y P )v y -rpv, V«GC°°(H), 

and using the fact that (compare ( 4.11[ )) 

p(Lw 1 )(x p ,y p ) = {LpWi){x,y). 

Prom (4.41), the Harnack inequality f?J Theorem 1.16] yields the estimate 

sup w\ < C inf ui\ < Cw(0), 
Bi/ 2 (o)n{j/ p >o} B 1/2 {o)n{y p >o} 

for a constant, C, depending at most on the coefficients of L. Combining the preceding inequality 
with (14.401) yields 



sup wi < Cw(Pq) < 20CNp, 



that is, 



sup wi < 20CNp, < p < p , 



(4.43) 



for a constant, po, depending at most on the coefficients of L. 

We will next bound w 2 on P+(Po), following the same reasoning as in the proof of Proposition 
Recall that < w 2 < w on B+(P ) by ( |4.39[ ) and w 2 = on in <9P+(P ) by Q4.37D - Assume 



4.1 



that Pi = (xi,y\) is a maximum point for the function w 2 on P+(Po) and that w 2 {P\) > 0. 
Therefore, Pi G P+ (P ) U r , where (by our convention) r = {y = 0} n dB+(P ). 

Case 1 (Pi G <f(u)). If Pi G S (u) n (P+(P ) U r ), then u(Pi) = ^(Pi). Recalling that 
w = u — ^ + £ by (4.35), we conclude from (4.34) and ( |4.39 ) that at Pi we have the bound 

w 2 {P x ) < w(Pi) = e(Pi) < 20iVp, 

provided p < pq. 

Case 2 (P x G "if (it)). If Pi G ^(n)n(P+(P )ur ) (where u > i/>), then since w 2 G ff 2 (S+(P ), »)n 
C(P p +(P )) obeys 

Pw 2 = a.e. on #(u) n S+(P ), 

by (4.4), ( |4.36 ), and ( 4.37| ), the regularity result in Proposition 3.5 implies that u>2 also belongs to 
C s 2+a (^(it) n (P+(P ) UT )). But w 2 achieves a positive maximum at Pi G 'jf(n) n (P+(P ) UT ) 
and so the strong maximum principle (Theorem 3.9) implies that w 2 must be constant on the 
connected component of ^(u) n (P+(Po) U To) containing Pi. Since W2 = 0onEn <9P+(Po) and 
w 2 {Pi) > 0, it follows that W 2 (P%) = w 2 (P 2 ), for some point P 2 with P 2 G <f (u) n (-B+(P ) U T ). 



4 See also [161 Theorem 4.5.3] for a version of the Harnack inequality for the linearization of the parabolic porous 
medium equation. 
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(Recall that, by hypothesis, Pq G &(u) and so ${u) n (P+(Pq) UTq) is non-empty.) We conclude 



that by fl4.34| ), ( |4.35p , ( |4.39| |, and the fact that u(P 2 ) = ^(^2), 

w 2 (Pi) = w 2 (P 2 ) < w{P 2 ) = £(P 2 ) < 2(Wp, 

provided p < pq. 

Combining the two cases and recalling that w 2 < w 2 (P\) on P+(Pq), by definition of Pi, yields 



sup w 2 < 20Np, < p < p . 



(4.44) 



Combining the estimates (4.43) and (4.44), respectively, for w\ and w 2 yields 



sup w < C'Np, < p < po, 



for a constant C which depends at most on the coefficients of L, and iV is as in (4.32). This 
yields the desired bound (4.30). □ 



Corollary 4.6 (Linear growth of solution near free and degenerate boundaries). Under the 
hypotheses of Proposition \4-4\ there are a constant < po < 1 and a constant < C < 00, 
depending at most on the coefficients of L, such that 



sup (u-V(Po)) < C P\\^\\c^{B+{p Q )y 0<p<p . 



(4.45) 



5. Proof of main theorem 

We will establish in this section the C^' 1 regularity of our solution, u, in Theorem For a 
much simpler example — interior C 1 ' 1 regularity for a solution, u, to min{Au — l,u} = on a 
bounded domain in M n — but one which conveys some of the flavor of our proof of Theorem 
see the proof of Theorem 1.1 in |19l p. 11], which is based in turn on ideas of Caffarelli [3]. 



Proof of Theorem \l.l\ Because of Proposition 3.10 we may assume without loss of generality 
that u G H 2 (V, to) D C(V) is a solution to the homogeneous obstacle problem (3.7) with obstacle 



function V G C 2+a {V) and / = on V. Recall that V = B+ o (Q ) is as in ([L5]), for some R > 
and Qq = (po,qo) G S. We may also assume without loss of generality that 



and also that 



< P < 1, 



A = 1 and < q < 1, 



since L is uniformly elliptic on V = Br {Qq) when qo > 1 and standard results imply that 
u G C^ X {V) [231 Theorem 4.38]. 
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Figure 5.1. Regions in the proof of Theorem 1.1 for estimating the Cl' norm of a solution. 
Let Pq = (xq, yo) G c to{u) n U, where U = B~^ o ^ 2 (Qo) as in (1.5) with R = Rq/2; see Figure 



Assuming without loss of generality thaljj 



5.1 



< y < 1, 



we will establish the bound 



y \D 2 u(P )\ + \Du(P )\ + \u(P )\ <C[\\u 



\C(V) 



+ 



(5.1) 



(5.2) 



where the constant C = C(L,Rq) may depend Rq and the coefficients of L. Since the constant 



C will not depend on yo (if Vo obeys (5.1)) this will provide the desired Ci' 1 bound on u up to 
y = 0. Set 



Vi := B. 



3i?o/4 



(Qc 



Since u is continuous on V, the exercise region, (o(u), as defined in (4.2), is a relatively closed 
subset of V U Tq . We may suppose without loss of generality that 



Otherwise, V\ C ^(u), where ^{u) is the continuation region, as defined in (4.1), and because 
Lu = on ^(u), Theorem |1 . 1| would follow immediately from Proposition 3.5 Now let d be the 
maximum number such that 

B dyo (P o )ng(u)nv 1 = 0. (5.3) 

Then there exists at least one poinl|^] 

Pi = ( Xl , yi ) G dB dyo (P ) n &(u) n Vi. (5.4) 



^Our assumptions so far that Po G U and Ro = 1 ensure < yo < 3/2, but standard results apply when yo > 1. 
®We alert the reader that in section [4] we use Po to denote a point in ,^F{u) whereas in this section we use Po to 
denote a point in ^(u) and Pi to denote a point in ,^(u). 



22 



P. DASKALOPOULOS AND P. FEEHAN 



Since P G U and Pi € Vi, we have < dy < 5Rq/A. 

Throughout this section, we let < po < 1 denote the smaller of the two constants in Propo- 
sitions 4.1 and 4.4 and, by replacing po with a smaller constant if needed, we may assume that 

< po < min{l,i? /5}. (5.5) 

We shall distinguish between three situations. We begin with the first situation. 




Figure 5.2. Regions for Step [I] (where d < po/A) of the proof of Theorem for 
estimating the Cl' 1 norm of a solution near the free boundary and away from the 
degenerate boundary. 



Step 1 {d < po/A). Since po < 1, we have B poyi (Pi) <g H (unless y% = 0), while Pi S V% implies 
dist(Pi,H n dV) > Po/4; since we also have po < Ro/5 by (5.5), we may conclude that 

5 



Payi 



[Pi)^V if < yi < 



(5.6) 



Since Pi = (xi,yi) G dB dyo (P Q ) and P = (x ,y ), we have \yi - y \ < y d < y /4 and thus 

3y /4 < yi < 5y /4, (5.7) 



and hence yo and y\ are comparable. In particular, we have < y\ < 5/4 by (5.1 ), and hence by 
(5.6) (see Figure 5.2), we see that 

B Poyi (Pi)mV. 



Set p := Ad and let ( G C 2+a (B pyi (Pi)) be the function defined by (4.9) (with Po and yo replaced 
by Pi and yi, respectively), that is, 



L( = L(l Pl ) onP m (Pi), 
( = WMy lP 2 on dBpy^Pi), 



(5.8 
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where (compare ( |4,8| )) 

M := m\\ cl , HBpQyiiPl)) , (5.9) 

and K is a constant which depends at most on the coefficients of L, and the following inequality 
holds (compare (4.7)) 



\L(\ = \L(l Pl )\ < M on B m (P x ) 



(5.10) 



It follows from (4.10) that 



My lP 2 <(< UMy lP 2 on B m {P x ). (5.11) 

Moreover, since Pi = (x\,y\) G dB dyo (Po) and 3yo/4 < y\ and d = p/A < 1/4, we have 
dist(Pi,Po) = dyo and 2dyo < 8dyi/3 = 2pyi/3 < py±, and thus (see Figure 5.2) 

B dyo (P ) C B^Px). (5.12) 



Therefore, C is also defined on B dyo (Po) and satisfies the bounds ( 5.11 ) on B dyo (Pq) with p replaced 
by d = p/4: 

^M Vl d 2 < C < ^M Vl d 2 on B dyo (P ), 



3 70 

— My d 2 <(< ^My d 2 on P rfyo (P )- 



and thus, applying (5.7), 



As in the proof of Proposition |4.1[ we set 

w ~ u - l Pl + C G i? 2 (£ m (Pi)) n 
The inequality (4.29) (with the role of B pyo / 2 (Po) there replaced by B pyi / 2 (Pi)) yields 

< w < CMy lP 2 on P m/2 (Pi), 



(5.13) 
(5.14) 



and thus, since yi < 5y /4 by ( |577| ) and p = 4c? and -6^/2(^0) C B pyi/2 (Pi) by ( |5.12[ ), 

0<^<C7My d 2 onP dW2 (P ), (5.15) 

for a larger constant C depending at most on the coefficients of L and where M is as in ( |5.9[ ) 
(compare (4.8)). 

Because B dyo (Po) C ^(w), we have Lu = on B dyo (Po), while L(l Pl ) = LQ on B pyi (Pi) by 



(5.8). It follows that 



Liu = onfl^fPo), 



(5.16) 



since B dyo (Po) C B pyi (Pi) by (5.12). Consider now the rescaled solution, w G C 2+a (Pd), given 
by 

w(x,y) :=w(x + y x,y + y y), V(x,y) E B d , (5.17) 
to the uniformly elliptic equation 

L yo w = on B d , (5.18) 
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where B d = Bd(0,0) and L yo is given by (4.12). The classical Schauder interior estimates for 
strictly elliptic equations pjfl Corollary 6.3] yield 

d\\Dw 



C(B d/4 ) + d 2 \\D 2 w\\ c[ B d/i) < C||tD|| c( B j, 



for a constant C depending at most on the coefficients of L (noting that d < 1/4). Combining 
the preceding inequality with the inequalities ( |5.15 ) for w implies the bounds 



\Dw(0)\ + \D z w(0)\ < -p\H\ C (B d/2 ) < CMy 



where C depends at most on the coefficients of L, and M is as in (5.9). Hence, since Dw(0) 
y Dw(P ) and D 2 w{0) = y 2 D 2 w(P ) by fl5.17[ ), we obtain 

^\Dw(P )\ +Vo\D 2 w(Po)\ < CMy . 

We conclude that 

\Dw(P )\ + yo\D 2 w{P )\ < CM, 
Similarly, the rescaled function ( G C 2+a (B d ) given by 

C(x, y) ■■= C(x + y x, y + yoy), V(x, y) G B d 



(5.19) 
(5.20) 



satisfies the uniformly elliptic equation (see (4.11)) 

L yoC = VoL( = VoL(l Pl ) = y /i on B d , 
where f\ := L{lp ± ) is a smooth, linear function with 

\\h\\c H B d )<c(mPo)\+m(p )\) 



< c\\n c ^(B dyQ {p )) 

= CM, 



< C 



C^(B P0V1 {Pi)) 



while C depends at most on the coefficients of L, and M is as in (5.9). Define rj G C 2+a (B\) by 
C(x,y) =: d 2 i](x/d,y/d), (x,y) := (x/d,y/d) G Bi, 

The function 7] obeys 



(L yo C){x,y) 



[Cxx + 2QCr(xy + cr 2 Cyy) V) + (r-q 



yo (i + y) 



Cx(x,y) 



and 



+ « (0 - 2/0 (i + 2/)) Cj/Oe) y) - r yoC(x, y) 

= ~Y^~ {Vxs + 2Qvmv + V) + d ( r ~ 1 ~ ~ ^ — ) ^(x, y) 

+ dK(6 - y (l + dy)) r/ 5 (x, y) - ry d 2 r](x, y) 
= : L yo4 rj(x,y), V(x,y) G S ls 

L yo ,dV{x,y) = yofi(x,y) = y fi{dx,dy) =: y fi(x,y), V(x,y) G B x . 
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We have Df\{x,y) = d~ l Df\(x,y) and so, noting that < d < 1/4, 

ll/l||ci(Bi) = ll/i|lc(Bi) + HAfillcCBi) 

= \\h\\c{B d ) + d\\Dh\\ c{Sd) <CM. 

Applying the classical Schauder interior estimates [121 Corollary 6.3] to the solution 7] to Lyg ^rj = 
fi on Bi gives 



\\ D v\\c(B 1/2 ) + \\ D v\\c(B 1/2 ) < C [ \\v\\c(B 1 ) +yo||/i||ci(Bi) 



for a constant C depending at most on the coefficients of L (recall that < uq < 1). Therefore, 
on B d , 

d-'WDCWciB^) + \\D 2 C\\c (Bd/2 ) < C (d- 2 ||Clc(B d ) + My Q ) , 
for a constant C depending at most on the coefficients of L. Combining the preceding inequality 
with the bound (5.13) for £ yields 

d^DCm + \D 2 C(0)\ < C (V 2 ||Clc(B d ) + My 

= c(d- 2 \\(\\ c{Bdyo{Po)) + My 
< CMy , 

for a larger constant C, but depending at most on the coefficients of L. Hence, since D(,(0) 
y D((P ) and D 2 ((0) = y 2 D 2 ((P ) by ( pO) ), we obtain 

d^yolDdPo) | +^|£> 2 CCPo)| < CMy , 
and thus, noting that 4 < d , 

|^C(Po)| + W)|^ 3 C(i%)l<CM. 



(5.21) 



Recalling that u; = u — l Pl + £ by (5.14), we conclude from (5.19) and (5.21) that 

\Du(Po)\ + ya\D 2 u(P )\ < CM, 



where M is as in (5.9) and so (5.2) holds for this step. 

We consider the second situation. 
Step 2 (d > 1). We shall consider two cases. 



Case 1 (d > 1 and dyo > po/4). Since dyo > /?o/4 for this case (see Figure 5.3) 

B + PQ „(Po) c B+ o (Po) nvc V(u) n v, 

and so Lu = on B^ o ^(Pq). The Schauder estimate (3.5) therefore yields 



Hlcy ( B+ /8 (P )) ^ C W U \\c^{B+ po /8 (P )) ^ ^Nl 



C(B p +/4( P 0))' 



(5.22) 
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Figure 5.3. Regions for Case [T] of Step [2] (where d > 1 and dyo > /5o/4) of the 
proof of Theorem 1.1 for estimating the Ci' 1 norm of a solution near the free 



boundary and near the degenerate boundary. 



with a constant, C, depending at most on a,po, and the coefficients of L, recalling that we have 
chosen < yo < 1 by our assumption (5.1) for this section (and thus A = 1). This yields the 



desired bound (5.2) for this case. 























} t 




d yo \ \ \ 
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Figure 5.4. Regions for Case [2] of Step [2] (where d > 1 and dy < Po/4) of the 
proof of Theorem 



1.1 



for estimating the C.}' 1 norm of a solution near the free 
boundary and near the degenerate boundary. (The radii 2dyo and po are not 
drawn to scale, since 2dyo < po/2 for this case.) 
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Case 2 (d > 1 and dyo < Po/4). Since cfo/o < Po/4, and also Pi E V\ and dist(Pi,Po) = dyo by 
(5.4) and po < Ro/2 by (5.5), we have (see Figure 5.4) 

^i ( p o) c B^jp,) c B+ (p 1 ) ^ur„, 



Thus, it follows from (4.45) (with Po replaced by Pi and p = 4dyo < Po) and Taylor's theorem 
(since dist(Pi,Po) = dyo) that 

sup (u - i/)(P )) < sup (u - ip(Po)) 

< sup (n-V(Pi)) + |^(Pi)-V'( J Po)| 



<Cdy 



CM(B+(F )), 



for a constant, C, depending at most on the coefficients of L, and hence 

sup (u - 4>(P )) < Cdy \\i}\\ c i,i { yy 

We now consider the function 

w:=u- V^(Po) E C 2+a (B+ (P )) n C(P+ (P )), 



(5.23) 



(5.24) 



which satisfies the equation 

Lw = -LiP(P ) = riP(P ) onP+ o (P ), 
since Bd yo {Po) fl <f (u) = 0. By defining the rescaled function 

w(x, y) := -j—w(x + dy x, y + dy y), V(x, y) E D, 
dyo 

on P := Pi n {y > -1/d} = Bi n {yd > 0}, we see that 

L d w = rifj(P Q ) on D, 



(5.25) 



with (compare (4.11) and (4.12)) 



Law := yd (w xx + 2gaw xy + 



<J 2 w yy ) + (r 



'—) 



w x + re (0 - dy y d ) w y - rdy w, 



and y d ■= 1/d + y. The operator L d becomes degenerate at y d = or equivalently y = —1/d 
which explains why the domain of consideration in the new variables is the intersection, D. 



It follows from the bound (5.23) that 



< C 



\\w\\c(B) ^ ^imicLMv-), 

for a constant, C, depending at most on the coefficients of L. Denote Di/ 2 ■= P1/2 H {yd > 0}. 
Hence, combining the preceding inequality with the Schauder estimate (3.5) for the solution w 
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to the equation (5.25) and noting that ya = 1/d + y > 1/d on {y > 0}, yields the bound 



D 2 w{0,0)\ + \DW{0,0)\ < \\y d D 2 w\ 



C(D 1/2 ) 



+ \\Dw\ 



C(D 1/2 ) 



<C \\w 



\C{D) 



+ rw(P ) 



recalling that w(Pq) = dyow(0,0); here, C is a constant depending at most on the coefficients of 
L. Since D 2 w{0,0) = dy D 2 w(P ) and Dw{0,0) = Dw(P ), we obtain 

y \D 2 w(P )\ + \Dw(P )\ < C|M| c i,i ( v), 



and thus, by (5.24), 



y \D 2 u(P )\ + \Du(P )\ < CUW^yy 



for a possibly larger constant, C, but depending at most on the coefficients of L. This implies 
(5.2) for this case. 

We consider the third situation. 

Step 3 (/Oq/4 < d < 1). This is the simplest situation. As in Step[2j we consider two cases. 




Figure 5.5. Regions for Case [T] of Step [3] (where p /4 < d < 1 and dy > Po/4) 
of the proof of Theorem 1.1 for estimating the Ci' 1 norm of a solution near the 



free boundary and away from the degenerate boundary. 



Case 1 (/)o/4 < d < 1 and dyo > po/4:). When dyo > /?o/4 (see Figure 5.5), we have Lu 
B~^ o , 4 (Pq), the estimate (5.22) for u holds, and (5.2) follows in this case. 



in 
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(5.26) 



Case 2 (po/4 < d < 1 and dyo < Po/4). We now assume that dyo < po/A (see Figure 5.6). We 
consider 

w:=u- V>(Po) G C 2+Q (£ %0 (P )) n C(B dyo (P )), 
and the rescaled function, 

u)(x,y) := — u;(xo + yo^,yo + yoy), (z,y) e 
yo 



which satisfies (compare (4.11)) 



Lw = rip(P ) on B d , 



with (compare (4.12)) 



Lw := (w xa; + 2gaw xy + <7 2 w yy ) + [ r 

+ n(6 - y (l + y)) u) y - ry w)- 



yo(i + y) 



The operator L is strictly elliptic on B d with ellipticity constant bounded below by a positive 
constant depending at most on the coefficients of L. In addition, since dyo < po/4, the bound 
(5.23) applies (irrespective of whether d < 1 or d > 1) to give 

\w\ < CdyQ||V>|| c i,i(y) on B dyo (P ), 

and thus 

H < Cdll^H^i,!^ on Srf, 
for a constant, C, depending at most on the coefficients of L. Combining the preceding estimate 




Figure 5.6. Regions for Case [2] of Step [3] (where p /4 < d < 1 and dy < Po/4) 
of the proof of Theorem 



1.1 



for estimating the Cl' 1 norm of a solution near the 
free boundary and away from the degenerate boundary. (The radii dyo and po are 
not drawn to scale, since dyo < p/4 for this case.) 
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with the classical Schauder interior estimate [.12, Corollary 6.3] gives 

\D 2 w(0)\ + \Dw(0)\ < \\w\\ c2+a(Sd/2) 

<c(\\w\\c(B d ) + \^(Po) 

< CU\\ c i, H v } , 

again for a constant, C, depending at most on the coefficients of L (recall that po/4 < d < 1 in 
this case and that po depends at most on the coefficients of L). Hence, 

y \D 2 w(P )\ + \Dw(P )\ < C|^|| c i,i(v), 
since D 2 w(0,0) = y D 2 w(P Q ) and Dw(Q,Q) = Dw(P ). Thus, by fljTggfr , 

y \D 2 u(P )\ + \Du(P )\ < C\\n c ^(vy 
for a possibly larger constant, C, and ( |5.2[ ) follows in this case too. 

This completes the proof of Theorem □ 
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